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temperature 
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SOLUTION TO THE DIFFERENTIAL EQUATION FOR COMBINED 
RADIATIVE AND CONVECTIVE COOLING FOR A HEATED SPHERE 

INTRODUCTION 

The simulation of space processing (i.e. , low gravity) cooling of metal 
droplets in free fall leads to a differential equation of heat transfer that has 
apparently not been solved in closed form [l]. Consider a small heated sf^ere or 
droplet falling vertically through a large diameter tube in which it is subjected to 
simultaneous cooling by radiation to the tube walls and forced convection to a 
gaseous medium. It is required that the time to reach a specific (cooler) tem- 
perature be explicitly determined. This will be done by solving the differential 
equation for time as a function of temperature. 


SOLUTION 


Let T denote the temperature of the sphere at any time t. The energy 
balance relating the heat transfer from the droplet by radiation and convection 
is expressed by 


= -€Aa(T^ - T ^) - hA(T - T ) 
dt o o 


(la) 


or 


= -K (T^ - T ^) - H(T - T ) 
dt o o o 


(lb) 


The first term on the right expresses the Stefan -Boltzmann law of thermal 

radiation; the second is the forced convection term. It is assumed that the 

coefficients do not vary during the short time of fall involved in the drop tube, 

which, in effect, constrains the problem to the class: slow cooling. Thus, 

K = cArr/mc and H = hA/mc. In addition, the droplet is considered so 
o 

small that any thermal lag between its center and surface may be neglected. 



The differential equation may be written in integral form as 


^ J_ r dT 


K (T - T )(T^ + T + T + T ^ + H ) 
o o o o o o 


+ const. (2) 


where H = II/K . Consider the cubic equation which may be formed from the 
o o 

cubic polynomial in the denominator 


T^ + TT^+T^T+T^+H = 0 


( 3 ) 


The three roots of this equation are; 



with 


Cl 


ki + 



Co 


2 





( 3 ) 


2 



and 


H 

ki = 20 + 27 

o (6) 

kj = n/ (ki)^+ 32 


The single fraction of equation (2) may now be broken into four fractions 
using the quantities rj, r 2 » and r 3 ; 


t 


X / f r 

K I ' T - T ' 
o \ o 


B dT 
T-n 



C dT 
“ + 
T-rj 


f 


D dT 
TT3 


+ const. 


(7) 


wheretA, B, C, and D are constants that are functions of T , rj, r 2 , and r 3 . 
Integration of equation (7) yields 


t = \<Aln (T - T ) + B /n (T - ri) + C ^n (T - rj) 

iV I o 

o L 

+ D^n(T-r 3 )J + const. (h) 

Using standard techniques from the theory of rational fractions, the constants 
i/L, B, C, and D are evaluated: 


(T^-r,)(T^-r3)(T^-r3) 


1 

(ri - - r2)(r, - r3) 


B 


(1*2 -T^)(r2 -ri)(r2 - rs) 


(rs “ T^)(r3 - rj)(r3 - rj) 

Substituting the system, equation (9), into equation (8) yields 


( 9 ) 


1 

= KjT^-r,)(T^-r,)(T^-r,) 

^ (ri - r2)(r, - rs) ” ^ 


-/n (T - T ) 
o 


(T^-r.)(T^-r,) 
(rs-r,)(rj-r,) - rj) 


(T„-r.)(T^-rs) , , 

" (^r-~(r~rs) 


+ const. 


( 10 ) 


Substituting the roots, equation (4), for the coefficients in equation (lO) results 
in 


(T^-ri)(T^-r2){T^-r3) = 4 (t/ + 


( 11 ) 


(rj - r2)(ri - r^) 


(Ct)3 - 4Ci -H 22 
3[(Ci)3 +2] 


( 12 ) 


(Tq ' ~ ^3^ ^ (C, + 4)[C,(Ci - 2) - 2NiTi(Ci + 1)1 

( ^2 - rj) ( rj - rj) iCjUCi)^ -f- 2J 


‘1 



( 14 ) 


(Tq - ri)(Tp - r;) ^ (c, + 4)[ Co(Ci - 2) + 2N/Ti(Ci + l)] 

(rg - ri) (rg - rj) ^ 


The roots r 2 and rg are complex conjugates, as are the coefficents of equations 
(13) and (l4). P rom complex number theory, the following formula will be 
used to evaluate the last two terms of equation (lO) which have conjugate 
coefficients and arguments, respectively: 

(u + Ivjin (z + iy) + (u - iv) fn (z - iy) = u in (z^ + y^) - 2v tan”^ ^ 

( 15 ) 


The conjugate roots are used to form 


T -T2 


T-rg 



(16) 


( 17) 


Using equations (13) through (17) , the following identifications of real and 
imaginary components are made; 


u 


3[(Ci)2+2] 


-2(C«4-l)(Ci + 4) 
■/3 C2[(Ci)2 + 2] 



-T C, 


\- = 


2 -JT 


(IH) 




Since the identifications, equation (18) , and equation (15) fit the form of the 
last two terms of equation (10) , proper substitution of the results of equation 
(ll) through (18) yields for equation (lO): 


t = 


K [4(T )^+ H 1 


- in (T - T 


o o 


. / [{C,)^-4C,*22] \ 

)' \3[{C^r*2]j 


in T 


+ -f (Ci + 1) 


(T^)'(C2)^ 

12 




-1 


V^^3■c^l(c,)’' + 21 




T C, 
o ‘ 


\ 


+ const. 


(19) 


This is the desired solution for the cooling time required to reach a 

specific temperature, T. The constant is evaluated in the usual way by using 

the initial temperature T. at t- 0, where, in general, T > T > T . 

i > i o 

For the case of radiative cooling only, the forced convection constant 

H = 0, so that C( = 2 and C 2 = ±2 ''/T, and equation (19) yields 


t 


1 

4K T ^ 
o o 



T + T 
o 

T - T 

o 


+ 


2 tan 



o 


+ const. 


(20) 


This is the same solution that the integration of equation (lb) gives with li = 0, 
thus affording a check of the solution, equation (19), for a simpler well-knowTi 
case. However, if in the solution, equation (19), we allow T 0, then the 
time for cooling is given by a limiting process as 



Cl 


(21) 


This same result is obtained if at the outset in equation (lb) we bet = 0, 

forming Uie Bernoulli [2] differential equation for this case whose solution also 
yields equation (21). This affords another check of the solution. 

CONCLUSIONS 


The solution to this differential equation will be applied as a guideline 
In the physics and engineering analysis of the cooling of small metal droplets 
;^!i.er solidification from molten metal as they fall freely down a tube. Such a 
fre«> fall will simulate for a few seconds the process of forming such droplets 
in a space laboratory. There are perhaps other physical processes in which the 
solution to such a differential equation could be useful. 
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